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We describe the resonant interaction of an atom with a strongly focused light beam by expanding 
the field in multipole waves. For a classical field, or when the field is described by a coherent state, 
we find that both intensity pattern and photon statistics of the scattered light are fully determined 
by a small set of parameters. One crucial parameter is the overlap of the field with the appropriate 
dipole wave corresponding to the relevant dipole transition in the atom. We calculate this overlap 
for a particular set of strongly focused longitudinally polarized light beams, whose spot size is only 
0.1A^ as discussed in S. Quabis et al, Appl. Phys. B 72, 109 (2001). 



I. INTRODUCTION 

For many applications, both in classical and quantum 
optics, it is a good idea to focus light to a small spot 
size. For example, in recent experiments on single atoms 
trapped in tiny dipole traps 0] the focusing system plays 
a major role. A natural question in that context is, how 
strong can the interaction between a single atom in free 
space and a light wave be? One expects the strength 
of the interaction to increase with decreasing focal spot 
size. However, there must be an upper limit to the in- 
teraction strength for at least two reasons: first, there 
is a limit to how much one can focus light of a given 
wavelength A; second, the cross section of a two-level 
atom is cr = 3A^/(27r), which seems to indicate that as 
long as all light is focused to within a the interaction is 
optimal. The latter argument, however, does not paint 
the complete picture. First of all, it leaves out polariza- 
tion effects, and second it does not distinguish between 
a classical object with a cross section a (such as a classi- 
cal oscillating dipole 0) and a quantum object with the 
same cross section (such asaJ = OtoJ=l transition in 
an atom). Here we take a closer theoretical look at the 
resonant interaction of a two-level atom with strongly 
focused light. 

To start with earlier theoretical work, in Refs. 0, Q 
a standard quantum-optical version of input-output the- 
ory (see 0j, and also j^) was used to calculate the ef- 
fects of light scattering off of a single atom in free space. 
The description of the light waves in the model has a 
1-dimensional character but one may expect a full 3- 
dimensional calculation to be necessary if the incoming 
light is strongly focused. Here we show how the input- 
output theory can be rephrased to include the full 3- 
dimensional description of the input and output light 
beams and yet keep all characteristics of the simpler 1- 
dimensional model. 

In other work 0, on the same subject exact 3- 
dimensional solutions of the Maxwell equations were con- 
structed by expanding the field in a complete set of func- 
tions that are well-suited to describe a cylindrically sym- 
metric beam. In the present paper we use a different 
method, and apply the results of Refs. 0, 0|, where the 



standard Debye approximation is used to construct ex- 
act expressions for strongly focused light beams. The 
authors find that a radially polarized beam of light pro- 
duces a field that is longitudinally polarized in the focal 
spot and that is focused down to a very small spot size, 
namely A « O.IA^. It is worth noting that this area is 
smaller than a by almost a factor of 5. 

One may wonder what the strongest focusing possible 
is. It turns out that that depends on one's definition. The 
smallest spot size is one criterion, but another, which is 
the most relevant for our discussion, is to find the max- 
imum possible electric field intensity in the origin given 
a fixed power for an incoming beam. That maximum is 
known to be achieved by an electric dipole wave (1L | . For 
illumination with a finite numerical aperture one still ob- 
tains the maximum by an electric dipole wave, but the 
optimum then depends on the polarization, see |il2j and 
also Section V. In any case, a measure of how focused a 
given light beam is, is thus given by its overlap with the 
appropriate dipole wave. 

Furthermore, for the interaction of light with an atom 
in the usual long- wavelength approximation |l3j| , one may 
expand the interaction in a multipole series. To lowest 
order, the atom interacts through an electric dipole inter- 
action; hence, to lowest order the only types of waves in- 
teracting with the atom are, again, electric dipole waves. 
Thus, it makes sense for two reasons to expand the field 
in multipole waves around the origin (which is where the 
atom is assumed to be). In addition, we note here that a 
multipole expansion is also convenient for the calculation 
of the field distribution in the focal region 14]. 

This paper is organized as follows: In Section II we 
review the essential properties of multipole waves, and in 
Section HI we consider the interaction of an atom with a 
field expanded in multipole waves. In Section IV we cal- 
culate two characteristic quantities of the scattered light 
as functions of the location of the photodetector: one 
is the intensity pattern and one is the photon statistics. 
The former quantity is classical in some sense, as a quan- 
tum object never gives rise to an intensity pattern that 
cannot be also obtained from the Maxwell equations and 
a classical scatterer. The photon statistics, on the other 
hand, does provide a quantum signature of the scatter- 
ing process. In Section V we consider the (longitudinally 
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polarized) focused light beams discussed in ^3 ^.nd 
calculate their overlap with the appropriate dipole wave. 
We compare the results to similar known results on trans- 
versely polarized waves. We conclude by raising several 
open questions concerning various aspects of focused light 
that are not treated in the main text. 



II. MULTIPOLE WAVES 

When quantizing the electromagnetic (EM) field, we 
may choose any complete set of mode functions to expand 
the electric and magnetic fields in. The usual choice is 
to take plane waves, but, as mentioned above, for the 
description of the interaction of radiation with atoms in 
the dipole approximation multipole waves are a good 
alternative. 

Multipole waves are eigenfunctions of commuting Her- 
mitian operators such that the excitations of such modes 
(i.e., the photons) possess definite amounts of energy E 
[E — huj], total angular momentum [J^ — J{J+l)h^], 
angular momentum in the z direction [Jz — Mh] , and 
parity P [either (—1)"'+^ (written here as P = X) or 
(-1)-^ (denoted by P = Z)]. We thus may expand the 
electric and magnetic field operators as 

cB{r) = 1 duj^Af^^u^^{r)a^^ + H.c, (1) 

where we abbreviated the set of discrete quantum num- 
bers as (J, M, P) = :vim&Y.u ■=J27=iJ2m=-jEp=x- 
Furthermore ajj^^ and a^ji, are the creation and annihila- 
tion operators for the mode with the corresponding eigen- 
values Lo and v. The functions $ are normalized to 



/ 



df$:,(r) • Cv'(r) = {2TTfS{u; - u;')S,,, 



(2) 



and the same for Here we used the obvious abbrevia- 
tion S,y,y' = Sjj'Smm'Spp'- The normalization factor 



2eo(27r) 



1/2 



(3) 



is chosen such that the free-field Hamiltonian takes the 
familiar form 
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r\E^{f) + c^B^{r) 



(4) 



The expressions for the multipole waves are somewhat 
complicated (see ^3); but fortunately we will only need 
the values in the origin (which is where the atom is as- 
sumed to be), and in the far field, at a distance r from 



the atom with rtOA/c ^ 1 with loa the relevant atomic 
resonance frequency. In the origin, the only waves for 
which the electric field (with which the atom interacts 
to first approximation) is nonzero, are the electric dipole 
waves. They are the waves with quantum numbers J ~ 1 
and even parity (corresponding to P = X). The other 
quantum numbers lo are arbitrary, and M takes on one 
of the values M = ±1, 0. In r = the dipole waves take 
on the values 



$(0) 



LulMX 



r.3/2 



Stt' 



1/2 



(5) 



where we explicitly displayed the dipole quantum num- 
bers V = {J =l,M,P = X). The unit vectors 



Mo 



z; u±i = {-iy^ x)/^/2 



(6) 



are the standard unit circular vectors. In the far field, a 
dipole field ^ujiMx{r) reduces to [Tj] 

1/2 

I .1 /i.n / — /I / X. I 1 1/1 



.cos{kor — 7r/2) 



c 



[uM - {uM ■ f)f], 

(7) 

with fco = uj/c. 

For later use it is also convenient to have the Fourier 
transform of a dipole wave at our disposal. It factorizes 
into the product of a function of the length of k and a 
function of the unit vector k = k/k, 

^1 / 3 

^koiMx{k) = 1— \/ ■^'5(fc - ko)[uM - {uM ■ k)k]. (8) 

fco V OTT 

Since we will consider (quasi-)monochromatic waves the 
K-dependent part is the most relevant for our purposes: 



^m{k) = \ -^[uM " [uM ■ k)k], (9) 
" Stt 



which is normalized to 



d2K|$M(K)|' = 1. (10) 



III. ATOM-FIELD INTERACTION 

Assume we have a two-level atom with a fixed dipole 
moment d = dux, with uk with K = 0,±1 one of 
the unit circular vectors. Through the standard elec- 
tric dipole coupling. Hint — —d-E, one gets the following 
interaction Hamiltonian in the rotating- wave approxima- 
tion 



Hi, 



ih / dwAt(w)[6^CT — cr^buj], 



where 



1KX\ 



1/2 



(11) 



(12) 
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Here we defined new creation and annihilation operators 
blj and b^^ describing the electric dipole wave that is cou- 
pled to the atom, with quantum numbers J — 1,M — 
K,P = X and arbitrary cu. We denote the discrete quan- 
tum numbers hy v — uq. Eq. Hill) is the same interaction 
Hamiltonian as that studied in Refs.[3,Q, but its validity 
goes beyond a 1-dimensional model by virtue of having 
expanded the EM field in multipole waves, as was already 
mentioned before in Ref. @. 

We can now write down the Heisenberg equations of 
motion for atomic and field operators. The equation for 
b^jit) is easily solved to give 

6^(t) =e-''^(*-*«)6oH + «H J dt'e-*"(*-*'V-(i'), 

(13) 

where to is a time in the far past [tg < t) and bo^uj) — 
buj{to)- The equations for the mode operators for all re- 
maining modes (uj, u) contain just the free evolution, so 
that 



a^,.W =e-''^(*-*«)ao(^,^)• 



(14) 



The standard Markov approximation now consists of as- 
suming that ti{uj) is more or less constant in the relevant 
frequence range around the atomic frequency to a, and we 
approximate 



r ^2, ,3 1 


1/2 


■ r ■ 


1/2 






2^ 





(15) 



with r the spontaneous emission rate constant. (Thus, 
one finds the same expression for F as in a plane-wave 
expansion.) Substituting this into the equations for the 
atomic operators gives 

&^ = -iuJA(T^ - ^o-" + Vfazhnit) 

= -r(l + (7,)-2Vfa+6in(t)-2\/fa-feJjt). 

(16) 

(The last equation also follows from — 2(t+(T^ ~ 1-) 
Here we defined the input (justifying its name from the 
fact that this field drives the atom) operator 5] 



&in(i) 



1 



dwe-^"(*-*«)6o(w), 



which satisfies the commutation relation 

{b,At)Mn^^{t-t'\ 



(17) 



(18) 



Substituting Eqs. 113|1 and (|14|l into the expansion of the 
electric field yields the usual separation into source and 
free parts, E = -Efroo + -Esourco- In the far field 



ir t) 

-^source V ' 



where for convenience we will only display the positive- 
frequency parts of the fields from now on; the negative- 
frequency part is just the Hermitian conjugate of the 



positive-frequency part, = (£?'^+^)'l'. Causality is 

obeyed as the atomic operator must be evaluated at a 
retarded time t — r/c with r = the distance from the 
atom. 

The free field contains two terms, one corresponding 
to the relevant dipole mode, the other to all remaining 
modes. The former depends only on the operator bo^Lu), 
the latter on ao{uj, v), according to (again valid in the far 
field) 



dujbo{u)e'''^'^'~*"^ cos{kr - tt/2) 



flLO 



1/2 



/67r. i/2.'Uk — {uk ■ f)f 
_2eo(27r);^J * r 

j dwao(w, iy)e-*"(*-*'') cos(/er - j7r/2) 



huj 



2eo(27r) 



'^\G^Y/2^^yPjM{fl^ (20) 



where Pji\i{f) is a transverse vector (perpendicular to 
r), whose form depends only on the quantum numbers 
J and M , and which is a function of the unit vector f 
|13| . Again making use of a Markov approximation, we 
can rewrite this in terms of the input fields, evaluated at 
the earlier time t — r/c, 

, UK - {uK ■ f)f 



where 



fibin{t - r/c)- 
r 3huJA nl/2 



r 

PjM{f) 



(21) 



(22) 



" 167reoc 

We introduced here another set of input fields (although 
here these fields do not drive the atom in any way). 



with commutation relations 

[aUt).<l{t')]=5{t~t'). 
The total field becomes then 

#+)(r,t) ^ [%/fa-(i-r/c) + 5i„(t-r/c)] 
UK - {uK ■ r)f 



(23) 
(24) 



Y <n(*-^/c)M 



PjM{f) 



(25) 



In the first line one recognizes the standard expression 
Hi Eq. (2.22)] for the output field operator 



bout{t)=bin{t)+VT(J-{t). 

This operator, too, satisfies 

[6out(i),6Lt(i')]=^(t-0- 



(26) 
(27) 
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IV. PHOTON FLUX AND STATISTICS 

In the following wc suppose we measure the flux of the 
output field and its statistics for a fixed polarization e 
at some fixed position r = i? in the far field. This is 
not an essential assumption, and we could easily define 
quantities similar to the ones defined below for every po- 
larization component detected. 

It is convenient to introduce yet another input field 
operator 



(28) 



where the presence of the geometric factor 

E I^^A/(i?)-ef (29) 

ensures that the relation 

[«in(i),aL(i')]=^(i-0 



(30) 



holds. Defining a similar geometric factor for the dipole 
field &, 



V = UK ■ e — {uk ■ R){R ■ e), 



(31) 



the detection operator can then be compactly written in 
terms of the operator 

C{t) = Va,,,{t)+V[h,,{t)) + ^/fa-{t)], (32) 



with R= \R\: namely, we get 



R/c). 



(33) 



We should note here that the operators appearing in C do 
not all commute. In particular, 5in(i) does not commute 
with the atomic operators, but ain(t) does. 

We can define two quantities of interest: a photon fiux 
operator (with the dimension of a rate) 



F^{C\t)C{t)), 



(34) 



and the second-order intensity correlation function at 
time zero, (/^^^(O), 



G(^)(0) 
F2 



G(2)(0) = {C^^{t)C^{t)). 



(35) 



The ordering of the noncommuting operators 6in and the 
atomic operators matters here. We can make use of a 
theorem given in that states that in expressions such 

as those for F and g^'^^Q) we can place Eitmcc to the left 
of e\^J^ and ^io"uLc to the right of E^;^]. That IS, we can 
place all operators to the right of atomic operators 
and b\„ to the left. 



This way, we can easily calculate these quantities in 
special cases of interest. In all cases we assume the field 
illuminating the atom has a central frequency ujl with a 
bandwidth B sufficiently narrow so that B is smaller than 
other rates in the problem, S <C F and B <C c/R. Of- 
ten, we will be interested in the steady-state (in a frame 
rotating at the laser frequency ojl) solution. Taking ex- 
pectation values and moving to a frame rotating at w^, 
transforms the Heisenberg equations (|f 6() into the optical 
Bloch equations 

{a-) = (zA-^)(a-) + x/f(a,&in(t)) 

(a,) - -r(l + (a,)) - 2Vf (a+6i„(i)) - 2Vf (6t„(t)a-), 

(36) 

where the laser detuning from resonance is A = ~^a- 
For later use we note that in the steady state the second 
equation gives 

Vf [(6^.7-) + (a+6i„)] + V{a+a-) - (37) 



A. Coherent states 

We are interested in the case of illumination with a field 
in a coherent state. One reason is that this corresponds to 
a laser field, another is that the statistics of the incoming 
light is then Poissonian, so that possible non-classical 
statistics arise from the scattering process, not from the 
incoming light. 

Hence, we assume the input fields satisfy 



6i„(t)|/3) = /3exp(-it^it)|/3), 
ain{t)\ot) = aeyi'p{—iujLt)\a) . 



(38) 



Thus, the input field is described by just two complex 
amplitudes, one for the relevant dipole part, one for the 
rest. We also define r/, a dimensionless number relating 
the total amplitude of the free field at the detection point 
to the contribution of the dipole field, as 



V-qP ^Va + Vj3. 



(39) 



Below we will connect these quantities to the overlaps 
with dipole and other multipole waves. The optical Bloch 
equations (|36|l depend only on the amplitude (3 



= -F(I + (a,))-2Vf/3(a+)-2Vf/3*(a-). 



(40) 



The steady-state solution is 



r a- 



~2/?(l + ^(5) 

1 + 6-' + myv 
i + s^ + 8\(3\yr' 



(41) 
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where we defined the dimensionless detuning S = 2A/r. 
It is now straightforward to calculate the quantities F 
and g(2)(o). 

First, we compare our result with that of Ref. 0,13 for 
resonant excitation ((5 = 0). We find for the flux F 



F (X 



(l-2/|r;|)2 + 8|/3|Vr 



1 + : 



Vr 



(42) 



Compare this with the expression from Eq. (30), 



F = n 



(1- 275/7)' + 87^7s/7^ 

1 + 8njs/7^ 



(43) 



with TZ the total incident flux, 75 the spontaneous emis- 
sion rate into the solid angle subtended by the incident 
beam and 7 the total spontaneous emission rate (cor- 
responding to our r). Extreme focusing corresponds to 
7s = 7/2. We can then make the identifications 



7^ 

7s/7 



(44) 



This comparison, though, is not perfect. In our case t] 
is complex and can, in principle, take on any value, but 
< '^is/l ^ 2 (and with light coming from 1 direction 
only, one even has 2^s/l < 1)- 

The strongest effects on photon statistics occur in the 
weak driving limit |/3p F, — in the strong driving limit 
the atom saturates and the output field will display Pois- 
sonian statistics — and on resonance. We therefore con- 
sider the special case of weak on-resonance excitation, 
and obtain 



F = \V\'\f3\'\fj-2\ 



(45) 



This result of f/'^^^(0) is plotted in Figure 1. Remarkably, 
the plot is very similar to that obtained in Ref. |8i] for 
illumination with Gaussian beams and detection in the 
forward direction (see Fig. 8 there). In the latter case, 
the result is plotted as a function of the beam waist. The 
remarkable aspect is, though, that the Gaussian beams 
are in fact no longer solutions to the Maxwell equations 
as the focusing conditions are too strong for the parax- 
ial approximation to be valid. Here, on the other hand, 
the result is not an approximation, but is plotted as a 
function of 

For reference we also give the flux of the output dipole 
field at position R in the absence of the atom: Fq — 
\Vf}\^. This leads us to the following observations in spe- 
cial cases 

77 = 1 If one illuminates the atom with only the dipole 
field, so that rj = 1, one has the strongest possible 
focusing. All light is coupled to the atom, and we 
note that the correspondence relation H44|l shows 
that this case corresponds to 275/7 = 2 in terms 




FIG. 1: (7'^'(0), in the weak-driving limit and on resonance, 
as a function of for various values of the complex phase 4> 
oi Tj: (j) = n/{2n) for n — 1 . . .7. 



of Carmichael's parameters, indicating illumination 
with light that fills the full 47r solid angle. More- 
over, this case is the only one where rj and thereby 
the flux and photon statistics do not depend on the 
detection point. 

For T] = I the output flux equals the input flux: 
(^out^out) (^L^in), and F = Fq. One way to 
understand this is to note that while (bin) = P, we 
have for the source field \/T{a^) = —2/3. Thus the 
expectation value of the total field is — /3. However, 
this does not mean that a tt phase change is the 
only difference. 



In fact, the statistics of the light has been affected 
by the presence of the atom. The output photons 
are bunched as g^^^'iO) = 9, even though the input 
field dispayed Poissonian statistics. The explana- 
tion for this effect is similar to that given in 0, |3 
for strong bunching: the atom cannot absorb any 
photons if it is in the excited state, but can if it is 
in the ground state. The detection of a photon thus 
makes it more likely the atom would be found in the 
excited state, which in turn makes it more likely to 
detect a second photon, namely one emitted by the 
atom. 

When the free fields interfere destructively at the 
observation point (so that rj = 0), the field 
there arises solely from the atom's fluorescence. 
That light, as is well-known, is anti-bunched and 
(;*^^^(0) — (as we are in the weak driving limit). 
(The reason is simple, the detection operator is 
C = 'D\/Va~ and applying this operator twice 
yields zero.) The flux F \s F = 4Fo as a result 
of the source field being twice as strong as the in- 
put dipole field, \/V{a-) = ~2p. 

2 When r/ = 2 one finds the largest bunching effect, 
with 5(2) (0) 00. This occurs simply because 
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the total intensity there vanishes, F ^ 0, as the 
source field (of amplitude —2/3) destructively inter- 
feres with the free field (with dipole wave and the 
rest each contributing an amplitude P). 

77 = 4 When ij ~ 4 the total detected flux is equal to 
F = 4^0 with now the free fields contributing 4/3 
to the amplitude, the source field subtracting 2(3, 
as before. Interestingly, we again have completely 
anti-bunched light, g'-^^O) = (at least in the 
low-intensity limit). Here is why (in a quantum- 
trajectory picture) the light is anti-bunched in this 
case: The detection operator is effectively C ~ 
27(4/3 + \/T<7^), since the state of the free ra- 
diation field (in the detection point) is a coher- 
ent state with amplitude 4/3. The steady state 
of the atom in between photodetection events is 
I"*/") = Id) ~ 2/3/-\/r|e) (valid to first order order 
in the small parameter P/y/T). After the first 
detection of a photon, we collapse the state onto 
C|V') oc |5)-4/3/Vr|e) =: The probability rate 
of another photon detection is proportional to the 
norm of the wave function C\(j)) = — 16/3^/-\/r|e), 
which is only of order r(|/3p/r)^. In words, the two 
paths to produce a photon after the first photode- 
tection event (one from the laser field, the other 
from the atom) interfere almost completely destruc- 
tively. Hence, g^^"^ (0) — > in the weak driving limit. 

We note here that this type of anti-bunching is con- 
nected to a collapse of the atom to the excited state 
after the first photodetection, in contrast to the 
case of pure fluorescence, where the atom is col- 
lapsed into the ground state. 

?7 ^ 00 If the field not coupled to the atom is large, the 
photon statistics of the total field will be domi- 
nated by that field, leading to Poissonian light with 
g(^\0)^l. 



V. THE DEBYE APPROXIMATION AND 
DIPOLE WAVES 

The overlap of the incoming field with the appropriate 
dipole wave is clearly the most crucial quantity, as only 
that part interacts with the atom. For instance, if we 
denote that overlap by Od, then it easy to see that the 
following two relations, involving parameters used before, 
hold: 



\v~i\' 



\Od\ 



(46) 



One question relevant in practice is, what is the largest 
overlap possible for a given numerical aperture? In or- 
der to answer this question we revert to the Debye ap- 
proximation. In that approximation the field in a focus. 



resulting from a high numerical aperture, is expanded in 
plane waves (hence, the resulting expressions are conve- 
niently given as Fourier transforms). The approximation 
consists in taking into account only the geometric op- 
tics rays from the aperture, thus leaving out edge effects. 
For a recent discussion of hig h-ap erture beams and more 
background information, see |1Q|. 

In Refs 9, 10] this approximation is used to calcu- 
late the intensity profile of a particularly strongly focused 
type of waves that were also generated in an actual ex- 
periment. Those waves are longitudinally polarized in 
the focal spot. As before, in reciprocal space the wave 
function factorize into a delta function 5{k — fcp) and a 
part depending only on the unit vector k. We only need 
the latter part, which we denote by In spherical coor- 
dinates (a, 13) for k one gets 



1 



'AT 

— otherwise 



A{a)p{a, (3) for a < 



(47) 



with p{a, f3) a unit vector indicating the direction of the 
electric field. 



p(a,/3) 



cos a cos [3 
cos a sin /3 
sin a. 



(48) 



The previous expressions all follow directly from Eq. (9) 
in Ref. 10]. Here is related to the numerical aperture 
by NA = smO, and the normalization factor M is given 
by 



1-9 

M=2tt dasina|A(a)p. 
Jo 



(49) 



The factor A{a) is determined by the input field on the 
lens. For the class of light beams studied in [gllTol]. 



A{a) = sin a ■\/ 1 cosa| exp(— sin^ a). 



(50) 



with a = f/wQ the ratio between the focal length / of 
the lens and the waist wq of the incoming Gaussian beam. 
The factor ^/\ cos a I is typical for an aplanatic lens sys- 
tem. 

In order to determine the overlap we need the Fourier 
transform of the relevant dipole wave. Here we need K — 
and $01 since the polarization in the origin is directed 
along Uf) = z. From Eq. we read off 



$0(0, /3) 



■ sin ap(a, /3) 



normalized to 







d/3 / dasina|$o(a,/3)|^ 
/o 



(51) 



(52) 



The overlap between x ^iid $o is then 

Od = 2'K I da sin Qx • *o = 27r / Aasir? aA{a)/VJI. 
Jo Jo 

(53) 
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FIG. 2; Dipole wave content p = lOdP, for the longitudinally 
polarized waves of Eq. 1501 for several values of a = f/wo'- the 
three bottom curves correspond to a = 2,1,0, respectively. 
The top curve gives the maximum possible overlap jdnaxp- 

This overlap has been calculated numerically as a func- 
tion of 9, and the results are plotted in Fig. 2. The 
plot also shows the maximum possible overlap for given 
numerical aperture. That maximum is achieved when 
A(a) (X sin a. This is most easily seen by defining a 
scalar product 

($,#):= / dasinQ!$-^*, (54) 

and noting one has to optimize this scalar product over 
all $ for fixed ^E" — sin ap. The maximum possible overlap 
with a dipole wave is then 

|Omax|' = ^ + ^cos3 0-^cos0. (55) 

In particular, for 9 = 7r/2, which corresponds to the 
strongest possible focusing given that the incoming light 
comes from one direction, one gets jOmaxP = 1/2, while 
for the types of beams considered here the best one can do 
is to go to the limit a ^ when jOdP = 64/147 « 43.5%. 

The figure also shows that for light coming from one 
direction the light beams from Refs. 0,0] are very close 
to the optimum (for fixed polarization) in the limit of 
small a, but for 9 > 7r/2 the distance from the optimum 
suddenly increases. 

In most optical experiments, however, light beams in 
the focus are transversely polarized. The overlap of such 
beams with a transversely polarized dipole wave (in our 
notation, such a linearly polarized dipole wave would, 
of course, be a superposition of $1 and $-1), has been 
studied before in Refs. 0,0. It turns out (see Fig. 1 
in 10 and Fig. 3 below) that for opening angles 9 less 
than 7r/2 the overlap with the appropriate dipole waves 
tends to be larger for transverse than for longitudinal 
polarizations. This is a consequence of the latter radia- 
tion pattern having larger side lobes. This observation 
is true both for a truncated dipole wave and for more 
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FIG. 3; Dipole wave content p — \Od\'^ for various waves. 
The dashed curve gives the maximum possible overlap with 
a transversely polarized electric dipole wave, obtained for a 
truncated dipole wave. The dotted curve gives the same but 
then for longitudinal polarization (and is the same dotted 
curve as in Fig. 2). The solid curve corresponds to uniform 
illumination with a transverse wave (called "Sine" wave in 
12]: the result given there is valid for 6 < 7r/2). These curves 
are all obtained from Ref. [T^. 

realistic light beams. For 9 = 7r/2 the dipole content of 
truncated dipole waves is exactly 50% for both polariza- 
tions, and for values of 9 larger than 7r/2 longitudinally 
polarized waves have larger overlaps. This is plotted in 
Fig. 3. Let us finally compare the numbers for more real- 
istic beams: The transversely polarized wave that results 
from (nearly) uniform illumination of an aplanatic lens 
(called "Sine" wave in 12]) has a dipole wave content 
= 32/75 « 42.7% at 6* = 7r/2, which is in fact 
slightly smaller than the number quoted above for the 
beams described by Eq. l|5n|l in the limit a ^ 0. 



VI. DISCUSSION 

We have shown how to reformulate the quasi 1- 
dimensional theory of 0, Q to study photon statistics 
effects in resonant light scattering off an atom in the full 
3-dimensional setting, by expanding the fields in multi- 
pole waves 0]. We applied this formalism to the case 
of illuminating an atom with strongly focused light, and 
found that there is a single parameter, 77, that contains 
all the important information. This is true when one as- 
sumes one detects a single polarization component. Oth- 
erwise, for each polarization component that is detected 
there is a parameter like rj. 

The overlap of the incoming light beam with an electric 
dipole wave determines the strength of the atom-light in- 
teraction. One expects that smaller focal spot sizes tend 
to correspond to larger overlaps, as indeed the dipole 
wave is the only wave with a nonzero intensity in the 
origin. We found a confirmation of this suspicion in Fig- 
ure 2, where the overlaps of a particular class of focused 
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light beams with very small spot sizes are plotted. These 
beams turn out to have almost the maximum possible 
overlap, given a fixed numerical aperture (this is true 
within the Debye approximation), and given a longitudi- 
nal polarization. With this measure, the more standard 
case of transverse polarization leads in fact to even bet- 
ter focusing for opening angles 9 less than tt/2, although 
the spot sizes is actually larger 0, For 9 > tt/2, on 
the other hand, longitudinal polarization becomes better 
than transverse polarization. 

Finally, we note that focused light contains a great 
deal more structure than just a small focal spot size. In 
particular there are several different types of topological 
properties that are both robust (i.e., they do not dis- 
appear when boundary conditions are slightly changed) 
and generic (i.e., they occur under general, non-special 
circumstances). For a nice discussion see 0|. Some 
of those properties, such as phase and polarization sin- 
gularities, occur on length scales at or even below the 
wavelength of the light. Similarly, around a zero of the 
intensity the spectral density of the light can be singular 
as well UM An open question is what the quantum sig- 
natures are of such topological structures, and whether 
one can probe those with an atom or a quantum dot? 
An atom is most likely too small to notice changes of 



polarization or spectrum, even if the change takes place 
within a wavelength, but a quantum dot may just be suf- 
ficiently large. In that case one may wonder, is there a 
difference between a quantum dot seeing red-shifted and 
blue-shifted light everywhere or seeing red-shifted light 
in one location and blue-shifted light in another? Or 
similarly, is there a difference between seeing unpolarized 
light or seeing horizontal polarization in one location and 
vertical polarization in another? 
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